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Abstract 

Based on recent studies of the temperature dependence of the energy and 
specific heat of Uquid nuclear matter, a phase transition is suggested at a 
temperature ~ .8 MeV. We apply Landau Ginzburg theory to this transition 
and determine the behaviour of the energy and specific heat close to the 
critical temperature in the condensed phase. 

The existence of an energy gap in the spectrum of even-even nuclei due to paired states 
of either protons or neutrons [|l| similar to that described by Bardeen, Cooper and Schrieffer 
(BCS) for electrons in a superconductor ||^ has led to the suggestion that nuclear matter 
should also exist in a condensed phase for some range of temperatures 0. The proper- 
ties of this superfluid phase in both nuclear and neutron matter have been studied in the 
BCS approximation using a variety of phenomenological forces as well as more realis- 
tic interactions Remarkably, all calculations yield qualitatively similar results for ^Sq 
pairing, namely that neutron matter exists in a condensed phase for kp less than about 
1.3 — 1.5 fm~^. Recent calculations, using the Paris potential 0, by the Catania group 
have shown that only slight deviations occur in nuclear matter. Such modifications, which 
can be characterized by the use of a smaller nuclear effective mass in the case of nuclear 
matter, are known to give rise to a slight decrease in the gap, A. Although such calculations 
suggest that such a low temperature phase should exist in both nuclear as well as neutron 
matter this has not been taken into account in, for example, astrophysical calculations since 
it is thought that it may be masked by other instabilities pf. 

In field theoretic language BCS theory is considered as the spontaneous symmetry break- 
ing of phase symmetry. The condensed phase, e.g. the superconducting phase, is charac- 
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terized by an order parameter, (A), which is zero at the critical temperature, Tc. It has 
been estabhshed that the order parameter and the critical temperature fulfill the following 
approximate relationship 

^ - 1.76 (1) 

where Aq is the value of the energy gap at T = and here we have taken the Boltzmann 
constant /c^ = 1- In the normal phase the order parameter is zero. Interestingly enough 
the same relationship has been found to hold in the aforementioned calculations in nuclear 
matter |^. Furthermore, it has been pointed out that the same relation between Aq and Tc 
also describes spontaneous symmetry breaking of chiral symmetry in QCD if Tc is taken to 
be 2/7^ ||lO|, where /^r is the pion decay constant. In all of the aforementioned cases the order 



parameter is obtained from a gap-like equation with appropriate quasiparticle interactions. 
Recent studies of nuclear matter have suggested that the origin of collective states may 



ultimately be linked to symmetry rearrangement [jTT|. This leads to a BCS-like condensed 
phase, separated from the normal phase, which has an order parameter that goes to zero at 
the critical temperature. Calculations in finite nuclei at finite temperature suggest that this 
provides a reasonable description of the vanishing of the collective degrees of freedom [|12| . 



Recently it has been demonstrated that the low temperature behavior of the specific 
heat of symmetric nuclear matter could be obtained from a finite temperature extension of 
the semi-empirical mass formula ||13[. The temperature dependence of the coefficients in 



the semi-empirical mass formula |jT^ was determined by fitting to the canonical ensemble 
average of the excitation energy of over 300 nuclei for temperatures T < 4 MeV, using 
experimental information of the energy spectra of nuclei in the mass region 22 <A< 250. 
The volume term was then used to determine the temperature dependence of the energy per 
nucleon and specific heat of nuclear matter. This displayed some rather interesting aspects: 
A structure in both the energy and the specific heat was observed at temperatures between 
.5 and 1.3 MeV (the structure in the specific heat is of course more pronounced). Below 
this temperature the behaviour of the specific heat was quite different from that expected 
for a Fermi gas of free nucleons . This is not unexpected as the low lying energy spectra 
of most nuclei are predominantly collective in nature. Above 1.3 MeV, the specific heat 
was essentially linear in temperature as is the case for a Fermi gas, but with the somewhat 
surprising feature that the slope coefficient was considerably larger than that suggested by 
the Fermi gas bulk level density parameter, 

1 m* , . 

a, ^ - — , 2 
15 m 



where it was assumed m* = (0.7 — 1.2)m |[T^. We propose that the larger slope might in fact 
be quite reasonable, based upon comparison with the case of liquid ^He. It is well known 
that at low temperatures, normal liquid '^He may be treated as a Fermi gas of quasi-particles. 
However, due to the strength of the interactions in the liquid, constants such as the quasi- 
particle mass are not easily calculated, and rather are determined from experiment. For 
liquid ^He an effective mass of m* ~ 3m [|16[ is obtained by measurement of the specific 



heat. Following a similar procedure, noting the similarity between the free interactions of 
^He and of nucleons, one anticipates that the normal liquid of nuclear matter would exhibit 
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properties that are quite similar to those of normal liquid He. In particular, an effective 



mass of m* ~ 2m is required to fit the specific heat given in [13|. It is important to note 
that in the treatment of liquid nuclear matter, it is usually assumed that the liquid may be 
replaced by a gas of quasi-particles with the mass of the quasi-particles being equal to that 
of free nucleons. Then the effect of nucleon-nucleon interactions leads to a new effective 
mass of the quasi-particle to be ~ 0.7 — 0.8. It is to be stressed that the procedure followed 
in the case of normal liquid ^He, namely to use the experimental data on specific heat to 
deduce the effective mass, is more satisfying. 

In this paper, we propose that there is a second order phase transition in liquid nuclear 
matter with a critical temperature and an order parameter rj. We apply Landau Ginzburg 
theory to determine the thermodynamic properties of the condensed phase close to Tc, from 
information about the normal phase. We follow the procedure used for liquid helium to 
determine the normal phase, namely modelling the system as a Fermi gas of quasi-particles 
with an effective mass m* ~ 2m, determined from the specific heat. We find that the 
behaviour of the energy per nucleon and specific heat across the phase transition with 
Tc ~ .8 MeV to be consistent with that shown in [|T^. 

Landau and Ginzburg have provided a simple theory of phase transitions which ap- 
proximates the free energy in the region around Tc and is most useful in analyzing the 
thermodynamics in this region. In particular, using only knowledge about the uncondensed 
phase one is able to make predictions about quantities in the condensed phase, such as spe- 
cific heat, magnetic susceptibility and compressability. Moreover, Landau Ginzburg theory 
can be derived from microscopic calculations . 

Following the Landau Ginzburg formulation it is necessary first to determine an expres- 
sion for the free energy per nucleon /(T) in both phases. In what follows the subscript 
1 will refer to the lower temperature (condensed) phase, and 2 to the higher temperature 
(uncondensed) phase. For the uncondensed phase, we take a quadratic form for the energy 
per nucleon which follows from a low temperature Fermi gas approximation of a normal 
Fermi liquid, 

W2iT) =a2 + k2T^, (3) 

where 02 and ^2 are constants. From the relations for the specific heat per nucleon in terms 
of W and the entropy per nucleon s, 

dW ds 

we are able to deduce the entropy per nucleon in the uncondensed phase, 

S2{T) = C2 + 2k2T, (5) 

where C2 is an unknown integration constant which later cancels out of the calculation. 
From eqs. (|]) and (|^) the free energy per nucleon in the higher temperature phase is given 
by 

f2{T) = a2-C2T~k2T^. (6) 
To determine the free energy per nucleon in the condensed phase, we make use of the 



Landau expansion for the free energy per nucleon in terms of an order parameter 77, 
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which goes to zero at the transition to the uncondensed phase. This order parameter is 
related to the presence of pairing expected at lower temperatures and vanishes with the 
pairing gap A at some critical temperature T^. The free energy per nucleon expansion to 
order r]^ is 

MT,v) = f2 + Ar]' + Bv\ (7) 

Here A and B are functions of temperature. The order parameter is determined by requiring 
the condensed phase to be stable below Tc (i.e. /i should be minimized w.r.t. rj). This leads 

to 

/. = A - (8) 

Further, since A is of opposite sign in the condensed and uncondensed phases, whilst B is 
strictly positive [^], the lowest order expansion of A in T — Tc can be parametrized as 



A{T)=a{T-T,)2jB{T,). (9) 



Note especially that a > is an essential requirement following from the phase diagram [jTT . 
Substituting for A{T), the free energy per nucleon near Tc is given by 

/i(T) = (a2 - aT^) + {2aT, - C2)T - {h + a)T\ (10) 

where /2 is taken from eq. (^). 

From the free energy per nucleon given by eq. (|T0|), we can now determine the energy 
per nucleon in the condensed phase near Tc, 

Wi{T) = (a2 - aT^) + (a + k2)T^ (11) 
= ai + hT^. (12) 

Comparing this to the uncondensed phase (eq. (|^)) we note that the T dependence is also 
quadratic, but has a larger coefficient. Thus the specific heat is discontinuous across the 
phase transition, and is necessarily larger {ki > ^2) in the condensed phase. 

We now compare the structure of Wi{T) and W2{T) to what has been determined from 



the finite temperature extension of the semi-emperical mass formula |T3|. Before proceeding, 
it should be noted that the energy per nucleon in nuclear matter is obtained from the volume 
term of the binding energy for finite sized nuclei. It may be anticipated that any sharp 
features (e.g. kink in the energy per nucleon and discontinuity in specific heat Acy) will 
appear smoothed out. Thus whilst comparison is still possible at a qualitative level, it is 
difficult to obtain quantitative estimates for the critical temperature and the discontinuity 
in the specific heat. 



In the region (.5 — 1.3 MeV) the energy per nucleon from |]T3[ is observed to show a peak 
above the simple T^ behavior. This is in good agreement with what might be expected 
from a smoothed out downwards kink in W at Tc, which follows from eqs. (^ and (p!2|). 
Furthermore, the specific heat in shows a sharp drop in the region (.5 — 1.3 MeV) which 
agrees well with a smoothed out discontinuous drop (= {ki — fc2)Tc). It should be noted 
that the specific heat above 1.3 MeV is very nearly linear, supporting the use of a quadratic 



4 



temperature dependence of W2, and that the slope below .5 MeV is greater than that above 
1.3 MeV, which is in good agreement with our prediction that ki > 

If we treat the uncondensed phase as a Fermi gas of quasi-particles with effective mass 
m* ~ 2m, as suggested earlier based on the linear behaviour of the specific heat, we can 
estimate this jump in specific heat at the transition to a condensed phase, where there is 
pairing with an associated energy gap A. This is given by 

Acv ~ 2A3cv (13) 

where cy is the specific heat per nucleon in the uncondensed phase. The uncondensed 
phase parameters for the energy per nucleon given in eqn (|^) have been fitted in giving 



k2 = 1/6.7 MeV~^ and 02 = —16 MeV. Using this value for k2 and assuming ~ .8 MeV we 
find Acv ~ .6 MeV, which is consistent with the behaviour of the specifc heat per nucleon 
given in [^. 

At temperatures considerably higher than Tc, the energy per particle given by eq. @ 
will become positive. It is reasonable to identify this with a transition from a Fermi liquid to 
a Fermi gas, at temperature Tlg- Using eq. (13) with the fitted parameters 02 and k2 taken 



from |T3[, we estimate this transition temperature to be at Tlg ~ 10 MeV. This compares 
favourably with Tlg 15 - 20 MeV in |1H[ and Tlg ~ 5 MeV (finite nuclei) 

In summary, following suggestions of a pairing transition in nuclear matter, we have 
applied Landau Ginzburg theory to provide estimates for the thermodynamical properties 
across such a phase transition. For information on the uncondensed phase, we assume a 
quadratic temperature dependence for the energy per nucleon. We find that the behaviour 
of the energy per nucleon and specific heat across the phase transition with T^ ~ .8 MeV 
is consistent with that observed from the fit of finite nuclei. Following an analogy with 
liquid ^He, we suggest that the essentially linear temperature behaviour of the specific heat 



observed in [1T3 is indicative that nuclear matter in the uncondensed (normal Fermi liquid) 



phase may be considered as a Fermi gas of quasi-particles with an effective mass m* ~ 2m. 
Such observation appears quite reasonable. 

Lack of experimental data on nuclear matter at finite temperature makes further refine- 
ment of the model difficult. Experimental determination of the thermodynamic properties in 
heavy ion collisions would be extremely helpful for understanding the properties of nuclear 



matter at finite temperature. Clearly the energy per nucleon obtained in at tempera- 
tures above 1.3 MeV is much stiffer than that of a Fermi gas of free nucleons, which is often 
used in many astrophysical calculations ^ which in turn should affect neutrino production 
rates in stars. As this is the major cooling mechanism in these objects it would be interesting 
to see precisely how important this deviation is. 
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